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Abstract 

A theoretical interpretation is given for the observed long-distance correlations in potential 
fluctuations in TJ-II. The value of the correlation increases above the critical point of the transition 
for the emergence of the plasma edge shear flow layer. Mean (i.e. surface averaged, zero-frequency) 
sheared flows cannot account for the experimental results. A model consisting of four envelope 
equations for the fluctuation level, the mean flow shear, the zonal flow amplitude shear, and the 
averaged pressure gradient is proposed. It is shown that the presence of zonal flows is essential to 
reproduce the main features of the experimental observations. 
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I. INTRODUCTION 



Transport barrier formation is mostly caused by the emergence of a radial electric field 
shear [H [21 E]- This radial electric field may be induced by poloidal flows and/or a gradient 
in the pressure, apart from the direct particle losses. 

A simple model for barrier formation and transition to a high confinement regime that 
was solely based on the poloidal flow shear was proposed in Ref. jl]. In this model a mean 
sheared flow is amplified by the Reynolds stress O |6l [7] and turbulence is suppressed by 
shearing [S]. The combination between those two effects allows having two possible types 
of states. On the one hand, states with vanishing flow shear and high turbulence level (low 
confinement). On the other hand, above a critical threshold, states with non-zero flow shear 
and reduced turbulence fluctuations (improved confinement). The transition between these 
two types of states is a continuous bifurcation. 

In Ref. [9] the model was extended by incorporating the pressure gradient component of 
the radial electric field. This extended model shows the existence of two critical points, the 
first being the same as in the previous model. The second transition, happening at higher 
density and temperature, is a discontinuous transition to a zero fiuctuation state where the 
radial electric field is only due to the pressure gradient. Experiments have shown [10] that 
the L to H transition [TTj leads to a high confinement state with the radial electric field 
shear dominated by the pressure gradient. That is why the second transition in this model 
has been associated to the L to H transition. 

Recently and in experiments carried out in the TJ-II stellarator, the first transition (linked 
to the generation of the poloidal flow) has been identifled [iSl [lH US] with the emergence 
of the plasma edge shear flow layer [161 HZI ■ 

New experimental results report the existence of long-range potential correlations 
in the toroidal direction. Two probes are set toroidally separated and not in the same 
fleld line. In addition, consider the intersection point of the fleld line going through the 
flrst probe with the plane of the second probe. The distance between this intersection 
point and the second probe is larger than a poloidal correlation length of the high-fc turbu- 
lence. These correlations are observed during the transition leading to the formation of the 
plasma edge shear flow layer in TJ-II. The observed correlations correspond to (non-zero) 
frequencies below 30 kHz and thus they cannot be explained by mean (i.e. surface averaged. 
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zero-frequency) sheared flows. In the present work we aim to show that the experimental 
findings of Ref. [T^ can be understood in the framework of simple transition models if one 
appropriately incorporates the contribution of zonal flows [T^]. Here we use the term zonal 
flow in the sense of low frequency fluctuating flows with = and small but non-zero kg. 
A transition model including zonal flows was proposed in Ref. [20], which we slightly extend 
here in order to interpret the TJ-II results. The structure of the model equations is based 
on quasilinear calculations from pressure-gradient-driven turbulence pT]. In this paper, we 
use a phenomenological approach, determining the main parameters of the model from the 
experimental results. We will see that the model is able to capture the essential features of 
the experimental observations. 

In Section II we present the transition model incorporating the effect of zonal flows. In 
Section III a comparison with the experimental data is performed. Conclusions are collected 
in Section IV. 



II. THE TRANSITION MODEL 

The model presented in this section is an extension of the one used in Ref. [13] to discuss 
the emergence of the plasma edge shear flow layer. This is a model formulated at a radial 
point. The dynamical variables are the fluctuation level envelope £ := ((n/no)^)^''^, the 
mean flow shear V := (9r(Ve), the zonal flow amplitude shear, Vzf := driVezp), and (minus) 
the normalized average pressure gradient A/" := —adr{p)/ (p)(0). Here a is the minor radius 
of the torus, (■) stands for the average over angle coordinates, and r = corresponds to the 
magnetic axis. The equations of the model are 

^ = ^.Af'/'S - a.AT-'/'S' - a^SAT-'/'iV' + V|f), (la) 
dr 

d,Ar-'/'s'v + d2Ar-'^'vlpV-bv, (ib) 



dr 

dVzF Oi 



Af-^/'^S'VzF + dzU-^'^£^V - 6VzF, (Ic) 



dr 1 + ^^AA-iV2 

70 

^ = f - D£J\f. (Id) 
dr 

The structure of these equations is based on a quasilinear approximation of resistive pressure- 
gradient driven turbulence (the resistive interchange mode, due to bad magnetic field line 
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curvature, is assumed to be the basic instability at the edge of TJ-II). The hnear eigen- 
functions and the dependence of the hnear growth rates on M were computed in [21]. In 
particular, we have used a fluid approach to calculate the poloidal velocity shear and the 
sheared radial electric field. The reason is that the TJ-II plasma edge, r/a> 0.8, is in the 
collisional regime. We would like to point out that for the range of powers and densities in 
TJ-II considered here, neoclassical theory is only applicable to the inner region, r/a< 0.25. 
In addition, it has been shown [22J that in this regime the ambipolar radial electric field in 
TJ-II is small and shearless for r > 5 cm and that the electron root [23] is the only accessible 
root. Therefore, in the edge region the fluid formulation seems to be adequate for the studies 
to be carried out in this paper. 



The first term on the right-hand side of Eq. (la) corresponds to the linear instability 
generation of turbulence, the second to the non-linear saturation of the instability, and the 
two last terms to the supression of turbulence by sheared mean and zonal flows. The first 



and second terms on the right-hand side of Eq. (lb) represent the generation of mean sheared 
flow by Reynolds stress, and the third one is the collisional damping term (analogous to the 



last term on the right-hand side of Eq. (Ic)). The two first terms on the right-hand side of 



Eq. (Ic) give the generation of zonal flow by Reynolds stress; in the first term the factor 
(1 + «2A/'^^V^/7o)^^ corresponds to the effect of zonal flow supression by a mean sheared 



flow. Finally, in Eq. (Id), DS is the anomalous particle diffusivity, and V the normalized 
incremental particle flux, the control parameter of the model. Diamagnetic effects in the 
momentum balance equation have been neglected because we consider small values of F. 
In terms of dimensionless variables, 

t = 7or, E=—8, V=.[^V, V^zf = ./^Vzf, N = U, (2) 
7o V 7o V 7o 



the equations read: 



dE 
~dt 
dy 
d^ 

dVzF _ Qi 
dt ~ 1 + N-W^ 

dN 



N^/^E - N'^/^E^ - N-^/^E{V^ + r^p), (3a) 
a^N-^'^E^V + a2N-^'^Vl^V - hV, (3b) 
AT-^/^eVzf + a-iN-^l'^E^V - bVzF, (3c) 



dt 



F - DEN. (3d) 
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where oi = 7001/0;^, 02 = 02/02, 03 = 7003/0^, b = 6/70, T = r/70, and D = D/ai. 



The form of the equation for the time evolution of zonal flows, Eq. (3c), coincides with 
the one proposed in Ref. [20], except for the term proportional to 03, which is absent in 
the latter reference. It is worth commenting on the physical origin of that term. In the 
framework of the paradigm of shear flow generation by turbulence the Reynolds stress gives 
a non-zero contribution when the turbulent eddies are distorted by the presence of global 
shear flows. If only the mean flow is present the Reynolds stress gives the mean shear flow 
amplification term that we have discussed in the past [1]. When, in addition, zonal flows 
exist, the Reynolds stress gives two main contributions to the zonal flow equation. One 
comes from the coupling of m and —m + q components of the eigenfunctions distorted by 



the zonal flow (the flrst term on the rhs of Eq. (3c)). The other comes from a similar 
coupling but with the distortion induced by the mean flow (the second term on the rhs of 
Eq. ([3c))). Here m is large and corresponds to the turbulent component of the flow, whereas 
q is the wave-number of the zonal flow. A detailed calculation of those terms will be provided 
elsewhere. 

As will be shown below, there is a qualitative difference between 03 = and 03 7^ 0. If 
03 = the model exhibits a continuous transition between the state with V = and the 
state with V ^ 0. In addition, the stable flxed points are such that Vzf = 0. However, if 
03 7^ the transition is discontinuous and the stable, improved conflnement state has both 
V and Vzf non-vanishing. Since the toroidal correlations will be associated to the existence 
of stationary zonal flows, it seems that small but non-zero is required. Obviously, for 
very small 03 it is not possible to directly (that is, according to the continuity of the order 
parameter) distinguish between a continuous and a discontinuous transition. 

A. The toroidal correlation 

In this subsection we will try to express the correlation of the potential fluctuations at 
two toroidal positions separated by a toroidal angle 6 in terms of the variables of our model. 
The formula deflning the correlation is 



$(r, e, ip, t) - e, if, t)) ](^r,e,^ + S,t)- ($(r, e,ip + S, t)) 

I ■ (4) 

%{r,9,^,t)~mr,9,v,t))^ )( (<l>(r, ^, ^ + 5, t) - (<l>(r, ^ + 5, t))) ) 
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Assume that a separation of time scales exists, so that one can write 

$(r, e, t) - ($(r, e, v,t)) = $zF(r, 6, t) + <l(r, 6, ^, t), (5) 

where ^zpif^yOyt) is related to the zonal flow and ^{r,6,ip,t) to high frequency turbulent 
fluctuations. The high frequency fluctuations have short correlation length in the toroidal 
direction, except when the positions are aligned with the field lines. Here we assume that 
this is never the case. The zonal flows are characterized by a low frequency and not having a 
dependence on the toroidal angle. Now, take a field line on the magnetic surface labeled by 
r passing through (6*, ip) and {9s, ip + 6). We are assuming that 6 is such that \r{9s — 9)\ ^ le, 
where Ig is the poloidal correlation length of the high-fc turbulence. 
For us, (■) denotes an average over 9 and ip. Eq. (|4]) takes the form: 

^ ($ZF(r,g,t)^) ^ 1 .g. 

($ZF(r, 9, m + (I.(r, 9, ^, tf) i + |gj^ ' ^ 

where we used that the toroidal correlation of turbulent fluctuations, 

(<l(r,^,<^,t)<l(r,^,<^ + 5,t)), (7) 

is negligible for 5 large enough. 

The scenario suggested by the above considerations in order to interpret the experimental 
results of Ref. [TS] is clear. We expect that in a ramping experiment crossing the critical 
point, (^'zf)/(^^) be zero below the critical point and non-zero above it. From Eq. ^ we 
deduce that this makes the correlation function, /i, grow during the transition. 

Let us finally write (|6]) in terms of the variables of the present model. Since the model 
equations can be derived from quasilinear calculations of a pressure-gradient-driven turbu- 
lence model (which in particular is a fluid model), we assume that the density perturbation 
is the result of the convection of the equilibrium density by the flow V = — V$ x B/i?^. 
That is, 

kednQ/dx 

Denote by {•} the spectrum average. Using that ($^) oc and 7^ oc N"^/^ (see 

Ref. [21]), we infer that 
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Also, VzF oc $ZF- Therefore, 



Hence, 



/$2\ ^2 

^^^7^7^' ^>0- (10) 



This is the formula we were looking for. It gives the toroidal correlation of the electrostatic 
potential in terms of the variables of our model. In particular, it shows in a manifest way 
that the zonal flow is responsible for the appearance of toroidal correlations. 

B. Fixed points 

The fixed points of Eqs. ([3]) are the solutions of 

N^/^E - N~^'^E^ - N-^/^E{V^ + V|p) = 0, (12a) 
{a,N~^/^E^ + a^N-^'^Vi^ - h)V = 0, (12b) 

+ - ^ZF + a,N''l'^E^V = 0, (12c) 

DEN = r. (12d) 
A fixed point corresponding to a low confinement regime always exists: 

(i) Vo = V^zFo = 0, Eo = {T/Dy/^\ No = {V / Dfl^\ 

Define the critical flux, Fc := D{ai/h)^^^^^ . The fixed point (i) is stable if F < Fc and 
unstable if F > Fc. 

It is easy to see that when F > Fc, there is another fixed point: 

(ii) Vo = 0, VzfI = N - N-y^E = {aJhf'^\V/Df"> - {ai/hy^\V I Dfl^"^ , Eq = 
{h/aif/^\TlDf'\ No = {ailhfl^\TlDf/\ 

which is always unstable. 

The discussion of the fixed point with 7^ is more difficult and depends on the value 
of 03. If 03 = and F > Fc there is an additional fixed point 

(iii) V^zFo = 0, V^o' = ^ - N-^/^E = (ai/6)3/iO(F/D)3/5 - (ai/6)-^/20(F/D)3/io, Eo = 
{b/aif/''>{T/Df/\ No = {ai/hf'^\T/Df'\ 



which is stable^. Note that the transition is continuous at F = Tc- A plot of the bifurcation 
is given in Fig. [T] 

When 03 = 0, the dynamics never reaches an equihbrium solution with Vzf 7^ 0. Equiv- 
alently, Vzf can be non-zero only transiently. But this is a problem for reproducing the 
long-range correlations observed during the transition to improved-confinement regimes in 
TJ-II. The results of Ref. |18j were obtained in ramping experiments traversing the critical 
point. In Fig. [2] we show the numerical results from our model when one performs a flux 
ramp traversing the critical point, going from a low confinement state to an improved con- 
finement one. When V starts growing from zero, and during a short time, Vzf follows it 
and becomes non-zero. However, at a certain moment the inhibition of Vzf by V becomes 
noticeable as V increases and Vzf decreases to zero after the transition. The interval of time 
in which Vzf 7^ coincides with the interval in which 7^ 0. However, in the experimental 
data one can see that the correlation has a non-zero stationary value above the critical point. 

The situation is very different for > 0. The fixed points (i) and (ii) and their linear 
stability remain unchanged. However, there is no solution with Vzf = and V ^ 0. The 
third fixed point, corresponding to the high confinement regime and which we will call 
(iii)', has both V and Vzf non-zero. In Figs. |3] and |4] we show numerical calculations of 
the bifurcation for different values of 03. The most remarkable feature is that for 03 7^ 
the transition is not continuous anymore, but the stationary values of the variables jump 
at Fc. Of course, the magnitude of the jump decreases when 03 decreases. Regarding the 
problem of long-range correlations, it seems essential to have non-zero 03. As shown in Fig. [5] 
this allows to have rampings in which the correlations reach a non-zero stationary value. 
Although there is no simple analytical expression for the fixed point (iii)' for a general value 
of F, we can give a good approximation for the supercritical solution at F = Fc. Define 
e •= and assume 

E = E, + e^E2 + o{€''), N = N, + e^N2 + o{e''), (13a) 
V = eV, + o{e^), Vzf = eVzFi + o(e'). (13b) 

Ec = {b/aiY^^ and Nc = b/ai are the stationary values of E and for 03 = 0. Introducing 
this expansion in Eqs. ( [I2| ) and solving for the lowest order we get: 

^ At least for moderate values of F. 
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2/3 

E2 = [ - ) (14a) 



o . , X 5/6 W13^ -20 P 



ai 



b ^ 



N2 = -{-j E2 (14b) 
13 fb\-'/' lOai fb^'^'' 



^2 10 ai / 6 



1/2 

^Fl = -^-(-) ^2. (I4d) 

At this point, we must comment on an issue. In the experiments, that are carried out by 
means of density ramps, it is difficult to distinguish between a continuous and a discontinuous 
transition. We can see that by comparing the averaged flows in Figs. |2] and [5] A continuous 
transition may appear very sharp because the velocity shear grows exponentially in the 
initial phase. The sharpness depends on the value of this growth rate and on the noise level 
from where the velocity emerges. 



III. COMPARISON OF MODEL AND EXPERIMENT 

Experiments were carried out in the TJ-II stellarator in Electron Cyclotron Resonance 
Heated plasmas (Pecrh < 400 kW, Bt = IT, {R) = 1.5 m, (a) < 0.22 m, t(a) G [1.5, 1.9]). 
The plasma density was varied in the range [0.35 ■ 10^^, 1 ■ 10^^]m^^. Different edge plasma 
parameters were simultaneously characterized in two different toroidal positions approxi- 
mately 160° apart using two similar multi-Langmuir probes, installed on fast reciprocating 
drives (approximately 1 m/s) [Mj- For details on the probe arrangement see Ref. jl8|. Probe 
1 is located in a top window entering vertically through one of the 'corners' of its beam- 
shaped plasma and at ~ 35° (where <f is the toroidal angle in the TJ-II reference system). 
Probe 2 is installed in a bottom window at ip ^ 195° and enters into the plasma through 
a region with a higher density of flux surfaces (i.e. lower flux expansion) than Probe 1. It 
is important to note that the field line passing through one of the probes is approximately 
100° poloidally apart when reaching the toroidal position of the other probe that is more 
than 5 m away. Edge radial profiles of different plasma parameters have been measured 
simultaneously at the two separated toroidal locations with very good agreement. Proffies 
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were obtained in both shot to shot and single shot scenarios with the two probes and in 
different plasma configurations. 

There are five parameters in the model and the A parameter in the determination of 
the correlation; this is apart from the input function F. The two parameters b and D 
are directly related to the dissipation terms, viscosity and transport. We determine them 
from the expected values of those terms at the plasma edge. We take the flow damping 
rate to be ua = 1.7 ■ 10^ s~^ and the particle diffusivity Dp = 10^ cm^s~^. Therefore, we 
used b = D = 0.1. The oi parameter is determined from the criticality condition and 
the experimental measure of the density profile at the critical point. The measured density 
profile at about the critical density in TJ-II is such that Nc ~ 1. Therefore, we take ai = 0.1. 

In this model the long-range correlations are controlled essentially by the ratio as/A. 



Observe that using Eqs. (13) and (14) we can find an approximate value of the correlation 
at the critical point: 

-(-^^)"' 

To have a reasonable level of correlation we need a^/X about 1/3. With the present data it 
is not possible to distinguish between the two parameters. Therefore, just for convenience, we 
have taken 03 = 0.01 and A = 0.03. Finally, the parameter 02 has not a very visible impact 
on the comparison with the data and we have chosen 02 = 0.5. 

The input function required in modeling each discharge is the flux function F(t). Since 
we are not doing a detailed modeling of data, but only a description of the main features, 
we have parameterized the flux using only linear dependences in time. A typical example 
of how the data are described by the model is shown in Figs. [6] and [7], corresponding to 
discharge 18229. This is a case with a ramp up and down where the plasma crosses the 
critical point twice, once in the way up and another in the way down. Similar results have 
been obtained for 10 discharges of TJ-II using the same set of values for the parameters. 
As one can see in Fig. |6] the parametrization of the flux only gives the main features of the 
experimental flux. The parameters of this linear function are determined by getting a good 
description of the density function, which in this case is the ion saturation current. In Fig. [7] 
we have plotted the ion saturation current. Fig. [7^, the averaged flow velocity shear. Fig. [7)d, 
the ion saturation current fluctuation. Fig. and the toroidal correlation. Fig. [7]i. The 
agreement between the experimental data and the model description is quite satisfactory. 
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especially noting the extreme simplicity of the latter. 

IV. CONCLUSIONS AND FURTHER WORK 

Recently, long-distance toroidal correlations in the electrostatic potential fluctuations 
have been observed in TJ-II [TH]. The value of the correlations increases above the crit- 
ical point of the transition for the emergence of the plasma edge sheared flow layer. In 
a previous work the transition was interpreted in terms of a simple model [13] consisting 
of envelope equations for the level of fluctuations, the mean flow shear and the averaged 
pressure gradient. 

In the present paper we have shown that the phenomenon of long-distance correlations 
requires the extension of the model so that the effect of zonal flows is taken into account. 
With the addition of an equation for the zonal flow amplitude shear the model is able to 
capture the basic features of the experimental results. 

The structure of the model equations is based on quasilinear calculations of resistive 
pressure-gradient-driven turbulence, and we leave for future publications the detailed cal- 
culations which should allow to compute some parameters of the model. Herein, we have 
determined those parameters by fitting the experimental data. 
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Figure 2: Time-evolution of the variables of the model for T{t) = 0.05 + 10 H, ai = b = D = 0.1, 
02 = 0.5, and 03 = 0. Top: E (sohd) and iV (dashed). Middle: V. Bottom: Vzf (solid) 
and fi (dashed). The initial conditions are ^(0) = lO^^, Vzf{0) = IQ-'^, E{0) = {T{0)/Dy/^^, 

N{o) = (r(o)/z))6/i3_ 
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Figure 3: Stable stationary values of V. The values of the parameters are ai = b = D = 0.1, 
02 = 0.5, and 03 = 0.01 (solid), 03 = 0.005 (dashed), = 0.001 (long-short dash). 
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Figure 4: Stable stationary values of Vzf- The values of the parameters are ai = b = D = 0.1, 
02 = 0.5, and as = 0.01 (solid), as = 0.005 (dashed), as = 0.001 (long-short dash). 
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Figure 6: Flux function, r(t), used in modeling the experimental data. 
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Figure 7: Comparison between the model equations Q and the experimental data for the shot 
18229. 
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